program mynormal_lf
	version 8.1
	args lnf mu lnsigma
	quietly replace `lnf' = ln(normden($ML_y1,`mu',exp(`lnsigma')))
end
